Quantum theory of photonic crystal polaritons 
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We formulate a full quantum mechanical theory of the interaction between electromagnetic modes 
in photonic crystal slabs and quantum well excitons embedded in the photonic structure. We apply 
the formalism to a high index dielectric layer with a periodic patterning suspended in air. The 
strong coupling between electromagnetic modes lying above the cladding light line and exciton 
center of mass eigenfunctions manifests itself with the typical anticrossing behavior. The resulting 
band dispersion corresponds to the quasi-particles coming from the mixing of electromagnetic and 
material excitations, which we call photonic crystal polaritons. We compare the results obtained 
by using the quantum theory to variable angle reflectance spectra coming from a scattering matrix 
approach, and we find very good quantitative agreement. 

PACS numbers: 42.50.Ct, 42.70.Qs, 71.36.+C, 73.20.Mf, 78.20.Bh 
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Introduction. Since the pioneering works of more 
than fifteen years ago by E. Yablonovitch and S. John 0, 
a great deal of research in physics has been dedicated 
to the study of the optical properties of photonic crys- 
tals 0- In the last few years much attention has been 
devoted to photonic crystal (PC) structures embedded in 
planar dielectric waveguides, which are also called pho- 
tonic crystal slabs. These are systems with a periodic 
dielectric response in the plane of the waveguide, while 
a dielectric mismatch is used to confine the electromag- 
netic field in the vertical direction. 

One of the main issues of PC slabs is the existence of the 
light line problem: only photonic modes lying below the 
cladding light line are truly guided and with zero intrin- 
sic linewidth, while the modes lying above the light line 
are radiative. A method to calculate the photonic band 
dispersion and the intrinsic linewidth (complex energies) 
both below and above the light line has been recently 
proposed 0,0]. Until now, however, the physics of PCs 
has been considered mostly for what concerns the "op- 
tical" point of view, that is studying the propagation of 
light in periodic dielectric media. 

In this work we analyze the radiation-matter interac- 
tion in PC slabs, by considering the effects of the in- 
terplay between the electromagnetic field and semicon- 
ductor quantum well (QW) excitons. In particular, we 
show that polaritonic effects are present in PC slabs 
when the exciton-photon coupling is larger than the in- 
trinsic radiative linewidth of a photonic mode above the 
cladding light line. The band dispersion is greatly mod- 
ified in the vicinity of the excitonic resonance, leading 
to the formation of mixed states which are analogous to 
exciton-polaritons in bulk semiconductors [5( and micro- 
cavities [6-8] and which we call photonic crystal polari- 
tons (PCPs). Experimental results showing the forma- 
tion of PCPs in organic-based systems were previously 
reported in Ref. 0. 

Theory. Here we briefly describe the quantum the- 
ory of PCPs, starting with the second quantized total 



hamiltonian of the system, which is given by 
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In Eq. 0), the first term indicates the photonic band 
dispersion (real part of the complex eigenenergies), a^n 
(^kn) being the destruction (creation) operators of a pho- 
ton with wave vector k and band number n. These en- 
ergies are obtained by expanding the magnetic field in 
terms of the guided modes of the effective waveguide, that 
is the one with an average dielectric constant, and then 
by solving the Maxwell equation as a linear eigenvalue 
problem |3( ■ The imaginary part of the photonic modes 
is calculated by using a perturbative approach U If we 
restrict our considerations to the system schematically 
shown in the inset of Fig. la, the direction of periodicity 
is x and thus the wave vector is given by the component 
k x . The modes in a PC slab suspended in air can be 
classified as even (odd) with respect to specular reflec- 
tion through the plane xy, and even (odd) with respect 
to the plane of incidence (that is xz in the case con- 
sidered here). Throughout this paper we consider only 
modes which are spatially even with respect to xy plane, 
and odd with respect to the vertical mirror plane (even 
TE modes). These modes interact with transverse QW 
excitons. The excitonic problem, whose solution gives 
the energies Ml kv , is treated by solving the Schrodinger 
equation for the exciton center of mass envelope func- 
tion in a periodic piecewise constant potential having the 
same patterning as the photonic crystal structure. The 
in-plane wave vector and an integer v are good quantum 
numbers for the exciton wavefunctions, due to the peri- 
odic potential and to the spatial dispersion of the center 
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FIG. 1: (a) Real and (b) imaginary parts of the complex photonic energies of the even TE mode as a function of the wave vector 
in the first Brillouin zone, for the system shown in the inset, without excitonic resonance (no QW); (c) real and (d) imaginary 
parts for the same structure with a QW at the center of the dielectric core, with an excitonic resonance at ftfio = 1.4 eV. The 
dashed lines in (a) and (c) represent the light line in air. 



of mass; the corresponding destruction (creation) opera- 
tors are b^ v (b^)- The exciton-photon coupling matrix 
elements, Cu_ ni/ , are calculated in terms of the microscopic 
physical quantities as 

C k „,= f^^) 1/ ^v,( r )|^E k „(r,).r,|0) . 

(2) 

In Eq. (J2J, is the all-electron exciton wavefunc- 

tion, Ek„ is the electric field profile for the photonic mode 
at frequency Wkn in the PC slab, and the sum is over all 
the QW electrons. In particular, CknK is found to de- 
pend on the overlap between the exciton envelope func- 
tion and the transverse electric field (for what concerns 
TE modes). It is proportional to {f/S) 1 / 2 , where f/S 
is the oscillator strength per unit area which depends on 
the QW thickness The hamiltonian is diagonalized 
by using a generalized Hopficld transformation to expand 
new destruction (creation) operators i\ as a linear 
combination of a^n (&£„) and b\^ v (&£„), with the con- 
dition [Pk,-ff] = -Ek-Ac _10]. The new eigenenergies 



correspond to mixed excitations of radiation and matter, 
i. e. the PCPs. 

All the results presented in this paper refer to a model 
structure, namely a high index dielectric core suspended 
in air with a one dimensional patterning, and a QW 
placed at the center of it (see inset of Fig. la). The 
dielectric constant of the semiconductor-based core layer 
is set to the value e — 12, the dielectric constant of the 
air being simply e a i r = 1- The slab thickness and the air 
fraction are set to the values d/a — 0.2 and r/a = 0.3 
respectively, the lattice constant being a = 350 run. A 
QW of thickness Lqw — 8 nm is considered, and the os- 
cillator strength per unit area is set to the typical value 
f/S = 8.4 x 10 12 cm~ 2 . The intrinsic radiative exciton 
linewidth is assumed to be T = 0.1 meV. 

Results and Discussion. In Fig. la we show the 
photonic band dispersion in the energy range 1.1-1.7 eV. 
In Fig. lb we display the corresponding imaginary part 
as a function of the wave vector, which shows a maxi- 
mum at about 1.4 eV. The imaginary part goes to zero 
at k x = 0, that is at normal incidence, and at k x = 0.74; 
this last behavior is due to the crossing of the light line, 
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FIG. 2: (a) Calculated surface reflectance of a TE polarized plane wave incident on the structure of the inset of Fig. la. The 
excitonic resonance is set at hQo = 1-58 eV; the angle of incidence is varied from 6 = 0° to 8 — 20°, with a step of 1°. (b) Real 
and (c) imaginary parts of the complex energies of the photonic crystal polaritons, calculated with the quantum theory. The 
squared points in (b) are extracted from the calculated reflectance curves shown in (a), and the polariton splitting is indicated 
by an arrow. 



as shown in Fig. la, corresponding to the photonic mode 
becoming truly guided and stationary. If the photonic 
imaginary part is larger than the exciton-photon coupling 
matrix element (which is of the order of a few meV), a 
QW placed at the center of the PC slab does not produce 
any important change in the photonic band dispersion. 
Indeed, this result is shown in Figs, lc and Id, in which 
we display the real and imaginary parts of the complex 
eigenenergies coming from the diagonalization of Eq. 
with an excitonic resonance at hflo = 1.4 eV. In this weak 
coupling regime the photon and the exciton are almost 
uncoupled, as it can be seen f rom the crossing of the 
two dispersion relations in Fig. lc. The imaginary part 
of the exciton increases by an order of magnitude corre- 
spondingly to the crossing point, but this has a negligible 
effect on the photonic radiative lincwidth, which is still 
an order of magnitude larger than the excitonic one (see 
Fig. Id). 

In order to observe the strong coupling regime, the en- 
ergy of the excitonic resonance has to lie where the imag- 
inary part of the photonic mode is smaller than the 
coupling matrix element. From Fig. la and lb we see 
that for M7o = 1.58 eV the imaginary part of the cor- 
responding photonic band is about Im(huj) = 1CP 3 eV. 
In Fig. 2a we show the results of a variable angle re- 
flectance calculation, which is done by a scattering matrix 
approach [ll| . The dielectric function of the QW layer is 
frequency dependent, with a resonance at hilo — 1.58 eV. 
It is known that the sharp features appearing in the re- 
flectance spectrum correspond to the excitation of pho- 
tonic modes above the light line 0, 0], thus giving a 



point (k, uj) of the corresponding photonic band disper- 
sion; the wave vector component parallel to the surface 
is given by k — (u>/c)sin0. The (k, u>) points extracted 
from the results of Fig. 2a are compared to the dis- 
persion of PCPs in Fig. 2b, in which the squares are 
the scattering matrix results. The quantum theory is 
in excellent agreement with the classical approach. In 
particular, both theories confirm the anticrossing behav- 
ior of exciton and photon modes, which is a clear ef- 
fect of the strong coupling regime. With the parameters 
used in this work, the polariton splitting is as high as 
10 meV at the anticrossing point k x = 0.14, as shown in 
Fig. 2b. This splitting is found to be slightly larger than 
in semiconductor-based microcavities 00- I n Fig- 2c, 
finally, the imaginary part of the PCPs complex eigenen- 
ergies is shown. At k x = 0.14 the imaginary parts of 
the upper and lower polariton branches become equal to 
the same value Im(-E) = 10^ 3 eV, thereby indicating that 
mixed states of radiation and matter form in the PC slab. 
The dispersionless curve at ha(E) = 10~ 4 eV in Fig. 2c 
corresponds to the uncoupled excitonic modes. 

Conclusions. We have shown that the quantum the- 
ory of the interaction between electromagnetic modes 
in photonic crystal slabs and QW excitons can describe 
both the weak and strong coupling regimes. Moreover, 
we have found a very good agreement between the quan- 
tum theory and a classical approach based on a scattering 
matrix method. We thus conclude that new mixed states 
of radiation and matter excitations can be experimen- 
tally measured in semiconductor-based photonic crystal 
slabs by using variable angle reflectance or transmittance 
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techniques, provided that the experimental conditions re- 
quired for being in the strong coupling regime could be 
satisfied. 
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